Stability of Chiral States, Role of Intermolecular Interactions and 

Molecular Parity Violation 



Mohammad Bahramlf] and Afshin Shafieelll 
Department of Chemistry, Sharif University of Technology, P.O.Box 11365-9516, Tehran, Iran 

We study the problem of stability of chiral states, also known as problem of chirality, within the 
framework of a two-dimensional approximation of a symmetric double-well potential. We show how 
the symmetry breaking of the potential due to the molecular parity violation can stop the tunneling 
in a coherent way, accordingly stabilize the chiral states. Then, we use the quantum Brownian motion 
within a linear Lindblad-type equation to model how the intermolecular interactions make the 
tunneling incoherent, thus inducing a racemization by dephasing. Finally, we investigate the normal 
physical conditions where the dephasing racemization does not suppress the effects of the molecular 
parity violation, accordingly the molecular parity violation may be observed experimentally. 
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I. INTRODUCTION 



Chirality is a fundamental concept in different 
branches of science. The synthesis of chiral molecules in 
enantiopure forms, in particular for pharmaceutical uses, 
is a great challenge in the synthetic and structural chem- 
istry [TH3]. In biology, perhaps the most outstanding un- 
solved problem is the origin of bimolecular homochirality 
[U [5] . Chirality is also connected to the molecular parity 
violation where the weak interactions result in tiny en- 
ergy differences between two chiral enantiomers [UMlO|. 
The chiral behavior of a molecule can be described by 
a particle of mass M moving in a double-well potential 
V{q) with two minima at q = ±qo/2, separated by a bar- 
rier Vq, as shown in Fig.l; here q is a generalized inversion 
coordinate [TTlfT3] . In the absence of weak interactions, 
V(q) is symmetric and so the corresponding eigenstates 
of Hamiltonian have definite parity and are delocalized in 
space. However, the chiral states are parity asymmetric 
and thus cannot be stationary states. The existence of 
delocalized eigenstates may contradict the usual chemi- 
cal view that considers the molecules as objects with a 
well-defined spatial structure. There are also strong ev- 
idences suggesting the existence of stable chiral states, 
e.g., many chiral molecules show stable optical activity 
which means the chiral states are stable. Accordingly, 
although the delocalized states of some molecules clearly 
confirm some experimental results (e.g., inversion fre- 
quency of ammonia jll] ) , for some other molecules there 
should be a mechanism to stabilize the chiral states. This 
leads to the problem of the stability of optical activity, 
known also as Hund's paradox or the problem of chirality 
[HI HSJ. Indeed, in the literature two different versions of 
chirality problem have been formulated. One is the origi- 
nal concern first expressed by Hund himself [HI [15] : why 
some chiral molecules show a stable optical activity, as- 
sociated to a well-defined chiral state, in spite of the fact 
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that these states are not eigenstates of the Hamiltonian, 
in particular none of them corresponds to the ground 
state. In recent years, this initial concern seems to have 
been transformed into a stronger one [TBIEH] : why chiral 
molecules have a definite chirality. The new version of the 
question has been extended to all non-planar molecules, 
chiral or non-chiral, since their potential energy surface 
has at least two minima [TTJ [T^] . This has led to a very 
lively debate. 

For an isolated molecule such as ammonia, the hydro- 
gens beat between two minima of V(q) with the frequency 
uj x . Using standard WKB methods, io x can be calculated 
as uj x = Aujo(qa/ ^ h/MuJo) exp [-B Vq/ojq], where A and 
B are constants, of order unity, which depend on the 
precise mathematical form of V (q) [13] ■ Obviously, the 
slower the beating (i.e., tunneling through the barrier), 
the longer the life time of an isolated localized molecule, 
if the molecule was initially prepared in a handed form. 
Accordingly, once produced, the chiral states of molecules 
with high inversion barrier are stable for a remarkably 
long time, as explained by Hund [HI [15]. However, upon 
a quantitative test, Hund's explanation seems rather in- 
sufficient for the molecules with low inversion barrier. 
Nor can Hund's answer explain why the habitual states 
of a molecule are the left-handed or right-handed isomers 
and not its parity eigenstates. After Hund, several at- 
tempts have been made to explain the problem by either 
taking into the account a small parity-violating term in 
molecular Hamiltonian [TUl [HI HE] , or by adding non- 
linear and stochastic terms to the dynamics because of 
interactions with a surrounding environment |17H28j . 

The molecular parity-violating energies (although not 
yet observed experimentally) may have possibly impor- 
tant consequences for the question of bio-homochirality 
[51 [TO] . They can also play a promising role in the un- 
derstanding and testing of the fundamental symmetries 
of physics [SlfTO] \29\ [30] . The molecular parity- violating 
energies has been also considered as a possible way to re- 
solve the problem of optical activity [TUl El HE] ■ Here we 
focus on the connection of parity- violation and the prob- 
lem of optical activity. We investigate this issue both 
for an isolated molecule and also for the one under the 
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FIG. 1: A double- well potential provides a primitive but very 
constructive and convenient description of chiral molecules. 
The characteristics of the double-well potential, in particular 
the height of barrier (Vo) and the separation of two minima 
(qo), can be determined from the spectroscopic data. The 
chiral states are localized at each minima. 
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FIG. 2: The fraction of molecules in the first (thick line) and 
second (solid line) eigen-energies of Ammonia (as an exam- 
ple of a molecule with double-well potential) as a function of 
temperature (K). The partition function is obtained from the 
spectroscopic data of Ammonia in [12]. 



influence of environmental interactions. 



II. QUANTUM DESCRIPTION OF CHIRAL 
MOLECULES 

. For all non-planar molecules (chiral or achiral), the 
potential energy surface has at least two minima [TTJ [T^] 
in which the corresponding nuclei configuration states 
(i.e., chiral states) are localized. Thus, the chiral 
molecule, in its simplest form, can be described by a 
one-dimensional Hamiltonian with a double-well poten- 
tial V (q). In the absence of weak interactions, V (q) 
is symmetric, and the energy eigenstates are also parity 
eigenstates: 

HiP n (q) = E n V>„ (?) , iPn (q) = V>n (-<?) • (1) 

with n = 1,2, 3, .... Taking into account the typical en- 
ergy spectrum of H, at the ordinary temperatures and 
at thermodynamic equilibrium, the fraction of molecules 
in the two lowest states, ^1,2 (g), are very close to unity 
(see Fig. 2). However, ipi^(q) are delocalized between the 
two minima of V (q), thus they cannot represent the chi- 
ral states of the molecule. The chiral states, which are 
localized each at one minimum, are then coherent super- 
position of ^1,2 (q)- 
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The set {ipLt V'-R; VVi>3} provides the physically natural 
basis for the Hilbert space of our problem. Using this 
basis, H takes a block-diagonal form since its matrix el- 



ements H 



i\ H \ j) are: 
Hm = HiR = Hli = HiL = 0, 



Hij = Ei 



'ho 



(3) 



where i,j > 3. We denote the 2x2 square matrix of H 
in the subspace of {f/'Li 4>r} by H . The matrix Hq is 
"self-contained", in the sense that the time evolution of 



ipL,ipR depends only on H . In other words, because of 
Eq. ([3| , i?o is sufficient for describing the behavior of the 
chiral molecule, when only the first two eigenstates are 
thermodynamically relevant (i.e., when Vq 3> kgT with 
ks the Boltzmann constant and T the temperature). The 
matrix elements of Hq can be written as 



Hq = UJq I + uj.ct 



(4) 



with a; = (u) x , 0, ui z ), 
and: 



c = (<7 X , a y , a z ) the Pauli matrices, 
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where u> z is the measure of molecular parity violations. 
The largest numerically calculated value of ui z is of the 
order of 3 x Kr 2 Hz (e.g., for S 2 C1 2 ) pHU)]. Although 
parity- violation effects are very small (compare it with Vq 
which is in the range 10 10 — 10 15 Hz LI2[H]) nevertheless 
when u z 3> lo x they dramatically change the eigenstates 
of Hq, from non-localized states tpi^ to localized states 
i/jl,r, which can be considered as a partial solution of the 
problem of optical activity. 

However, a chiral molecule is not an isolated system 
and it interacts constantly with its environment, e.g., 
through the collision with the background gas. Mean- 
while, one of the main challenges in observing the molec- 
ular parity violation is how to get rid of intermolecular 
collisions [7] . Since the parity violating terms are so small 
(at best, uj z is of the order of few Hertz [7]) the question 
may arise that in what experimental conditions (pressure 
and temperature) the molecular parity violation still de- 
termine the dynamics of chiral molecule and so may be 
observed experimentally. Our approach is different from 
the usual view manifesting effects of intermolecular colli- 
sions in frequency shift and broadening of the molecular 
lines. In the following, we discuss our position in detail. 
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III. THE EFFECT OF INTERMOLECULAR 
INTERACTIONS 

. Several mechanisms have been proposed to describe 
how intermolecular interactions affect the dynamics of 
chiral molecules [TTHzTTl |2"5H2"5] . Here we focus on inter- 
actions due to collisions (also known as the collisional 
decoherence) , which also is very common in nature. In 
the two-state approximation, the collisional decoherence 
master equation is given by a Lindblad-type master equa- 
tion P2HIU ESI : 



dp 
dt 



-i[H ,p] - ^ (p- 



(7) 



where A is the decoherence rate. An important property 
of Eq. (|7|) is the dynamical transition of an initially delo- 
calizcd superposition state (e.g., V'1.2) into an incoherent 
mixture of ipL and ipn when A 3> 0J X , to z [T?H2"T1 HZrJ] , thus 
results to racemization, sometimes called as racemiza- 
tion by dephasing [18]. However, the dephasing racem- 
ization should be distinguished from the typical thermal 
racemization. In thermal racemization the molecules are 
energetic enough to surmount the barrier Vq, while the 
incoherent tunneling through the barrier induced by in- 
termolecular interactions results to the dephasing racem- 
ization. We should remind that the decoherence process 
does not induce the chiral selectivity, i.e., it cannot be 
used to separate a racemic mixture into its chiral compo- 
nents or it does not result the enantiomeric excess if we 
start from a racemic mixture. 

The great challenge is determination of A which de- 
pends on the details of the intermolecular interactions. 
For example, in [55] the collisional decoherence effect on 
D2S2 molecules by Helium atoms as the bath particles 
has been numerically computed where they considered 
the chirally sensitive EQ-ED/ED-ED London dispersion 
interaction between a chiral and an achiral molecule. 
Here we use a more general and simple method, first in- 
troduced in [31], to compute the collisional decoherence 
rate for non-planar molecules in gas medium. 

The influences of intermolecular interactions on dy- 
namics of chiral molecule can be treated by a kinetic 
description based on interaction events that can be de- 
scribed as scattering processes. Thus, we may apply the 
theory of quantum Brownian motion [32] to describe the 
effects of collisional interactions. An important observa- 
tion is that ^1,2 in a chiral molecule are spatial super- 
positions of an atom or a group of atoms between two 
wells of double-well potential. We may consider the chi- 
ral molecule as a Brownian particle in a superposition 
of distinct spatial positions separated by the distance go- 
The Brownian particle interacts with a bath of scattering 
particles. We work in dilute gas limit where three-body 
collisions or correlated scatterings are negligible. Accord- 
ingly, when the scattering is recoil-free and also isotropic, 
the decoherence rate can be obtained as follows 



A = r t 



P 

dpv(p) —a(p,q ) 
m 



(8) 



where T t is the total collision rate, m the mass of bath 
particles, n the density of bath particles, p the bath par- 
ticle momentum, v (p) the distribution of such momenta, 
/ (p, 9) the scattering amplitude, 9 the scattering angle, 
and 
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with O = (2pqo/H) sin (0/2). We assume that, for all 
important bath particles' momenta p, the a(p,qo) does 
not change appreciably. We can then write: 



A ~ r t - n — a(p,q Q ) 
m 



(10) 



with p = \Z2mksT the most probable momentum of an 
ideal gas. As custom, we can distinguish limiting cases of 
Eq. (10). In the high-temperature limit where pqo h, 




we have sin 0/0 
rate one obtains: 

A~r, 



. Then, for the decoherence 
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For g of f ne order of few Angstrom (i.e., molecular 
ranges), and m ps 4 — 30gr/mol (from Helium to air), 



one gets: 



10- 



10 K. Accordingly, in nor- 
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mal conditions, a good estimation can be obtained if we 
consider the mean rate of collisions as a measure for the 
order of magnitude of decoherence rate. 

The total collision rate is determined by the total 
scattering cross-section. By exploiting the optical theo- 
rem, the total scattering cross-section is connected to the 
imaginary part of scattering amplitude / (p, 9) [35 ] 136 ) . 
Using the phase-shift methods, / (p, 9) can be determined 
as a function of interaction potential [351 136j . A simple 
and general expression of the total collision rate can be 
obtained in the hard-sphere approximation [35H37] : 
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(12) 



where P is the pressure of bath particles, r\ is the effec- 
tive hard-sphere of bath particle and r 2 is the effective 
hard-sphere of chiral molecule. 

We now investigate the regime in which the molecular 
parity violation is greater than the dephasing racemiza- 
tion effects of an environment. To give an estimate, let us 
consider a chiral molecule in the air. Noting that at best 
w z ~ 3 x 10 _2 Hz for heavy chiral molecules where recoil- 
free assumption is valid (i.e., m chiral moloculc 3> m a ir)j r 2 
has the range of a few Angstrom 37J, r*i ~ cin 2 is the ef- 



fective hard-sphere of Nitrogen gas, we find from Eq. ( 12 ) 
that in the limits: 



P(Pa) 



< 3 x 10 
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the parity-violating interactions are physically impor- 
tant (i.e., u> z > A) and they control the dynamics of 
chiral molecules. It is clear that at the normal con- 
ditions (P = latm and T = 300K), the decoherence 
rate is much greater than the strongest molecular parity- 
violation (Ap = iatm, t=3ook ~ 10 9 — 10 10 Hz) and so the 
environmental effects smear out the effects of parity- 
violating interactions. Introducing the temperature of 
few Kelvins, we obtain the pressure should be less than 
~ 10~ 7 Pa. These conditions are in a very good agree- 
ment with the experimental conditions in supersonic 
beam molecular spectroscopy which is the main tech- 
nique that is being used, aiming at measuring the molec- 
ular parity violation [7]. 

IV. CONCLUDING REMARKS 

. We show that how parity-violating interaction can 
stabilize the optical activity of an isolated molecule. Ac- 
cordingly, for some molecules the parity violating in- 
teractions can partially resolved the problem of optical 
activity. However, the question remains open why the 
chiral molecules have not been observed in a superpo- 
sition of their chiral states. To find a proper explana- 
tion, we should also include the environmental interac- 
tions (in particular collisions with a background gas) into 
the dynamics. The environment may induce the dynami- 
cal destruction of the quantum coherence, thus results to 
the racemization by dephasing. The important issue is 



to find the experimental conditions in which the parity- 
violating interactions are physically important, e.g., they 
are greater than racemization effects of intcrmolecular 
interactions. Considering the quantum coherence in chi- 
ral molecules as a spatial superposition, we used the 
quantum Brownian motion to compute the conditions in 
which the parity violating effects may have physically ob- 
servable effects. However, here we used the simplest rep- 
resentation of the interaction potential (i.e., hard sphere). 
More sophisticated model potential can be used to cal- 
culate the scattering amplitude, and so to arrive at more 
accurate quantitative results. 

The above calculations may also shed light on the 
conundrum of the origin of life and particularly origin 
of homo-chirality. There are some suggestions that the 
molecular parity-violating energies (although not yet ob- 
served experimentally) may have possibly important con- 
sequences for the question of bio-homochirality, and thus 
the origin of life jTHHl [HI EH- However, much more the- 
oretical and experimental work is needed to find out 
whether or not there is any connection between parity 
violation and bio-homochirality because in these stud- 
ies the chiral molecules should be considered as an open 
quantum system interacting with its environment. 
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